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We study two types of models for the superconducting layered nitride /3-MNCl(M=Hf,Zr); a single 
band model on a triangular lattice, and a two band model on a honeycomb lattice. We find that the 
former model does not suffice as an effecitive model, while the latter one can be a good candidate. 
We propose from the study on the two band model a possibility of spin fluctuation mediated d + id'- 
wave superconductivity in the doped /3-MNC1. We show that the relatively high T c obtained in the 
doped band insulator is a characteristic feature of the spin fluctuation mediated superconductivity 
on a honeycomb lattice. We also find that the gap anisotropy on the Fermi surface strongly increases 
upon increasing the doping concentration, and the inter layer hopping suppresses superconductivity. 
These results are in qualitative agreement with the experimental findings. 
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I. INTRODUCTION 

Layered nitride /3-AfNCl 1 (M =Hf, Zr) doped with 
carriers is one of the most interesting group of supercon- 
ductors. The mother compound /3-MNC1 is composed of 
alternate stacking of honeycomb (HC) MN bilayer and 
CI bilayer^ This is a band insulator and becomes a su- 
perconductor upon doping electrons by Na or Li interca- 
lation. They exhibit relatively high T c up to ~ 25K for 
M =Hf and - 15K for M =Zr. The bilayer HC lattice 
structure consisting of M and N is considered to be play- 
ing the main role in the occurrence of superconductivity, 
and the two dimensional nature of the superconductivity 
has been revealed by nuclear magnetic resonance (NMR) 3 
and muon spin relaxation (/xSR) studies4^ 

Despite the relatively high T c , experimental as well 
as theoretical studies indicate extremely low density of 
states (DOS) at the Fermi level£~— In fact, they have the 
highest T c among materials with the specific heat coeffi- 
cient 7 as small as ~ lmJ/molK 2 . The electron phonon 
coupling is also estimated to be weak,^r— and the iso- 
tope effect is found to be smalL 11 ' 12 These experiments 
suggest that some kind of unconventional pairing mech- 
anism may be at work, but on the other hand, NMR 
knight shift measurement suggests spin singlet pairing, 13 
and the tunneling spectroscopy^ and specific healji ex- 
periments find a fully open, seemingly s-wave like gap. 
However, regarding this fully open gap, recent experi- 
ments show that the anisotropy of the gap increases with 
doping. Namely, when the doping concentration is small, 
the specific heat coefficient increases linearly as a func- 
tion of the magnetic field, suggesting an isotropic gap, 15 
while a steep increase of the coefficient at low magnetic 
field is observed for higher dopingii 5 - Also, 2A/(fcsT c ) 
determined from the specific heat measurement ranges 
from ~ 5 in the lightly doped regime to less than ~ 3 
in the heavily doped regime^ which may be an indica- 
tion that the maximum and minimum values of the gap 
have a substantial difference when heavily doped, and the 
specific heat is mainly governed by the minimum value. 
Moreover, a recent superfluid density measurement by 
/iSR shows that the gap is nearly isotropic when the dop- 



ing concentration is small, while the anisotropy increases 
for large dopingJ^ Very recently, absence of coherence 
peak in the spin lattice relaxation rate has been found in 
an NMR experiment^ again suggesting unconventional 
pairing, most probably with some kind of sign change in 
the superconducting gap. Furthermore, for LLjHfNCl, an 
intercalation of organic molecules tetrahydrofuran (THF) 
between the layers is found to enhance T C J£ 

Given these experimental circumstances, here we study 
two types of models for the doped /3-MNC1, a two band 
model on a HC lattice, and a single band model where 
the nitrogen site is effectively integrated out. While the 
latter model does not seem to suffice for explaining the 
relatively high T c in this series of material, the former 
one provides a possibility of spin fluctuation mediated 
d-wave superconductivity in the doped /3-MNClJ^ Spin 
fluctuation mediated pairing in a doped band insulator 
may sound surprising and odd at first sight since one 
cannot expect strong spin fluctuations. Nevertheless, we 
show that the relatively high T c despite the low DOS 
and the temperature independent spin susceptibility is 
a characteristic feature of the spin fluctuation mediated 
superconductivity of a doped band insulator on the HC 
lattice. The most probable pairing state below T c is the 
d + id! state, whose gap anisotropy on the Fermi surface 
is found to be strongly enhanced upon increasing the car- 
rier concentration. We also study the effect of the three 
dimensionality, where the inter layer hopping suppresses 
superconductivity. 



II. FORMULATION 

We first obtain the effective model for /3-MNC1. Al- 
though the material consists of bilayer HC lattice, here 
we consider models on a single layer lattice, and try to 
reproduce the bands that lie close to the Fermi level in 
the first principles band calculations. The simplest one is 
a single band model on a triangular lattice, in which the 
nitrogen site degrees of freedom is effectively integrated 
out, so that there is only one site per unit cell. This 
is a model that corresponds to the single band Hubbard 
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model for the cuprates, where the oxygen 2p x , 2p y de- 
grees of freedom is integrated out from the original three 
band d-p model. By considering hopping integrals up to 
fourth nearest neighbors (FigfTJ upper panel), we can re- 
produce the single band that crosses the Fermi level in 
the first principles band calculations . 8 ' 10 ! 20 ' 21 The band 
dispersion of the single band model with t = 0.65 (eV), 
t-2.lt = -0.03, t 3 /t = 0.06, t 4 /t = -0.025 is shown by the 
green dash-dotted line in FigfTJ 




FIG. 1: (color online) Upper : the single band model on a 
triangular lattice. Middle : the two band model on a hon- 
eycomb lattice. Lower : the first principles band calculation 
taken from ref.— and the band dispersion of the single band 
model (dash-dotted green) and the two band model(dashed 
red). 

To go a step further, we also consider a model on a 
single layer HC lattice consisting of alternating " M" and 



"N" orbitals with a level offset, which roughly reproduces 
the conduction and valence bands closest to the Fermi 
level, as shown by the dashed lines in the middle panel of 
FiglU Here we take the nearest neighbor (M-N) hopping 
t = 1.2eV, the level offset A/t — 2.7, and the next nearest 
neighbor (M-M) hopping t'/t = 0.35. 

The on-site interaction U/t = 6 is considered in the 
single band model, and the on-site interaction U/t = 6 is 
introduced on both M and N orbitals in the two orbital 
model 2 ^. The band filling n is defined as the number of 
electrons/number of sites. For the single band model, 
n = 2x, where x is the Li or Na content, since there 
are two M sites per unit cell, which results in bonding 
and antibonding bands, and the electrons are doped only 
in the bonding band. For the two band model, n = 1 
corresponds to the non-doped case, and x and n is re- 
lated by n = 1 + x. We use fluctuation exchange (FLEX) 
method^—, which is kind of a self-consistent random 
phase approximation, to obtain the Green's function. 
Then we solve the linearized Eliashberg equation, 

k' V ,m' 

x V lm (k - k')Giv (k')G mm , {-k')<j> Vrnl (k'). 

Here, G is the Green's function matrix (with l,m,etc... 
labeling sites in a unit cell) obtained by FLEX. V is 
the spin singlet pairing interaction matrix given by V = 
|^ 2 Xs P — \U 2 Xch with the spin and charge suscepti- 
bility matrices Xsp(ch)(<?) = Xin-(?)[1 - (+)Uxhr {q)]~ l , 
where Xi rr (?) is the irreducible susceptibility matrix 
Xirr(g) = -j? Efc G ( k + l)G(k) (N is the number of k- 
point meshes). In the following, the maximum eigenvalue 
of the spin susceptibility matrix will be referred to as 
the spin susceptibility. We take up to 128 x 128 fc-point 
meshes and up to 16384 Matsubara frequencies. The 
eigenvalue of the Eliashberg equation A increases upon 
lowering the temperature and reaches unity at T = T c . 



III. SINGLE BAND MODEL: CALCULATION 
RESULTS 

We first present calculation results for the single band 
model. In Figf21 the eigenvalue of the Eliashberg equa- 
tion is shown as functions of temperature for spin singlet 
and triplet pairings for n = 0.2. It can be seen that 
the triplet pairing strongly dominates over singlet pair- 
ing. The gap of the spin triplet pairing has an /-wave 
form, where the nodes of the gap do not intersect the 
Fermi surface (the ridge of the Green's function squared) 
as shown in Fig|3l On the other hand, the singlet pairing 
has a d-wave gap (not shown). The possibility of this 
type of spin triplet /-wave pairing was in fact proposed 
on a canonical triangular lattice in refj 2 ^. The triplet 
pairing dominates over singlet pairing because ferromag- 
netic spin fluctuations arise due to the dilute band filling 
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FIG. 2: (color online) Eigenvalue of the Eliashberg equation 
as functions of temperature for /-wave and d-wave pairings 
in the single band model. 




FIG. 3: (color online) The contour plots of the FLEX result 
at the lowest Matsubara frequency for the single band model 
in the hexagonal Brillouin zone with n = 0.2 and T = 0.02t 
(a) The Green's function of the upper band squared, (b) the 
spin susceptibility, (c) the superconducting gap. 



the pairing occurs in the spin singlet channel^, we be- 
lieve that the single band model does not suffice as an 
effective model for /3-MNG1. 

IV. TWO BAND MODEL : CALCULATION 
RESULTS 

A. Superconducting transition temperature 

We now move on to the two band model. In this model, 
we find that spin singlet pairing strongly dominates over 
triplet pairing. In FigJH we plot T c of the singlet pairing 
superconductivity as a function of the band filling. It can 
be seen that T c ~ 30K is obtained, which can be consid- 
ered as relatively high noting that T c ~ 100K (assuming 
t ~ 0.4eV appropriate for the cuprates) is obtained by 
the same method for the Hubbard model on the square 
latticed The difference between the single and two band 
models is dramatic in that not only the T c but even the 
leading pairing symmetry is different. This is in con- 
trast to the case of the cuprates, where the single band 
Hubbard and three band d-p models is expected to give 
roughly the same conclusions. 
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as shown in FigFJJb), and also because the nodes of the 
/-wave gap do not intersect the Fermi surface. 

Nonetheless, as can be seen from FigJU the tempera- 
ture at which the /-wave eigenvalue reaches unity, if any, 
seems to be very low. We have also tried a model that 
considers the nearest neighbor off-site interaction V in 
addition to U. In this case, strong charge fluctuations 
arise, which makes the triplet vs. singlet competition 
more subtle, but in any case, the eigenvalue of the Eliash- 
berg equation remained small even when the magnitude 
of V is close to the point where the charge susceptibility 
diverges. Considering the experimental fact that the T c 
is relatively high in the doped /3-MNC1, and also that 



FIG. 4: Tc plotted as a function of the band filling for the two 
band model on the HC lattice. 

The high T c obtained in this model can be understood 
as a combination of several features of this model. Let 
us start with a canonical HC lattice, where t' = and 
A = 0. We show in FigG^a) the FLEX result of the 
maximum value of the spin susceptibility as a function 
of temperature for the band filling of n — 1.08. The spin 
susceptibility is nearly independent of T , 15 i 27 which is in 
sharp contrast with the Hubbard model on the square lat- 
tice. For example, for the square lattice with n = 0.7, we 
have a strong enhancement of the spin susceptibility upon 
lowering the temperature. With further hole doping to 
n = 0.65, the spin susceptibility is suppressed, but even 
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FIG. 5: (color online) (a) The maximum value of the spin sus- 
ceptibility (blue, left axis) and the eigenvalue of the Eliash- 
berg equation (red, right) as functions of temperature for the 
square lattice with n — 0.7 (dashed) or n — 0.65 (dash dotted) 
and for the honeycomb lattice (solid) for n = 1.08. U — 6t 
in all cases. HC and SQ stand for the honeycomb (with 
t' — 0) and the square lattices, respectively. The tempera- 
ture at which A = 1 is the T c . (b) T c as functions of the 
level offset A obtained by FLEX+Eliashberg equation for the 
square lattice with U — 6t, n — 0.7 (dashed) or for the HC 
lattice with U = 6t and n = 1.08 (solid). 



there, the spin susceptibility moderately increases upon 
lowering the temperature. Also in FigJSJa) , we show the 
eigenvalue A of the linearized Eliashberg equation as func- 
tions of temperature for the HC and square lattices. T c is 
the temperature where A reaches unity. The DOS at Ep 
is nearly the same for the square lattice with n — 0.65 
and the HC lattice with n = 1.08, and also the spin 
susceptibility has similar values at low temperature, but 
still, the HC lattice has higher T c 22. 

So what is the origin of this high T c ? FiglBJa) shows 
the contour plot of the Green's function squared, whose 
ridge corresponds to the Fermi surface, which consists 
of disconnected 2D pockets centered around K and K' 
points. The spin susceptibility in FigJJJfb) is maximized 
at wave vectors that bridge the opposite sides of each 
pieces of the Fermi surface. As can be seen in Fig|7J{c) 
(and more clearly in (d)) the gap has a d-wave form, 
where the gap changes sign across the wave vector at 
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FIG. 6: (color online) The contour plots of the FLEX result 
at the lowest Matsubara frequency for the honeycomb lattice 
(t' = A = 0) in the hexagonal Brillouin zone with n = 1.08 
and T — O.Olt (a) The Green's function of the upper band 
squared, (b) the spin susceptibility, (c) the superconducting 
gap of the upper band, (d) the Fermi surface (the two circles) 
and the sign of the gap function are schematically shown in 
the extended zone scheme. The dashed arrows represent the 
wave vectors at which the spin fluctuations develop. 



which the spin susceptibility is maximized. In contrast 
to the case of the square lattice, one of the nodes of the 
d-wave gap does not intersect the Fermi surface because 
of its disconnectivity, i.e., the gap is like "p-wave" if we 
focus only on one of the Fermi surfaces. Since smaller 
number of gap nodes on the Fermi surface is favorable 
for superconductivity^ this can be a reason for high T c 
despite the low DOS and the weak spin fluctuations. 

We now introduce the M-N level offset A as in the 
model for /3-MNC1. As shown in FigJSJb), we find that 
superconductivity is relatively robust against the intro- 
duction of A. This again is in sharp contrast with the 
square lattice, where finite A strongly suppresses super- 
conductivity. The weak effect of A may be because the 
DOS of the HC lattice is small near the band center, so 
that the introduction of A, which opens up a gap at the 
center of the band, is less effective than in the case of 
the square lattice, where the DOS diverges at the band 
center for A = 0. 

Let us now comment on the origin of the difference be- 
tween the single and two band models, namely, low T c 
triplet pairing in the former model and the high T c sin- 
glet pairing in the latter. The difference mainly comes 
from the difference in the band filling, i.e., for the single 
band model, the band filling is dilute, favoring ferromag- 
netic spin fluctuations, while in the two band model, it 
is near half filling in favor of antiferromagnetic spin flue- 




tuations. Of course, if the level offset A is large enough, 
the two band case should tend to the single band case, 
as have been shown in ref.—, but in the present case, the 
magnitude of A is small enough to retain the nature of 
the A = situation. 



B. Gap function 

Apart from the T c , another important issue is the form 
of the superconducting gap. By symmetry, there are two 
degenerate d-wave gaps, and the most probable form of 
the gap below T c is the form d + id! , where the two d- 
wave gaps mix with a phase shift of tt/2.—~— Since the 
two d-wave gaps have nodal lines at different positions, 
this kind of mixture leads to a gap that has a finite abso- 
lute value on the entire Fermi surface. The d + id' form 
of the gap can be constructed from the d-wave gap ob- 
tained by solving the linearized Eliashberg equation. In 
the upper panels of Fig 13 we show the contour plot of the 
d + id' gap for the case of n = 1.06 and n = 1.16 for the 
model for /3-MNC1. To clearly see the anisotropy on the 
Fermi surface, we plot the magnitude of the gap (normal- 
ized by the maximum value) in the lower panel of FigJT] 
It can be seen that the d + id' gap is nearly isotropic 
for low doping concentration, while it becomes more and 
more anisotropic for higher doping. The d + id' pairing 
scenario is consistent with the experimental findings that 
the superconductivity occurs in the spin singlelji 3 -, s-wave 
like channel,— ' 14 ' 15 while the coherence peak in the spin 
lattice relaxation rate is absent; 17 It is also consistent 
with the experimental finding that the gap anisotropy 
increases upon doping i 15 ' 16 Since the d + id' state breaks 
time reversal symmetry, it is interesting to experimen- 
tally investigate such a possibility. 



C. Effect of dimensionality 

Finally, we consider the effect of dimensionality by 
adopting a three dimensional model where the two band 
model on the HC lattice is connected by vertical (M-M 
and N-N) hopping integrals t z (inset of Fig[8]). Here, we 
take up to 32 x 32 x 32 fc-point meshes and up to 8192 
Matsubara frequencies. In Fig0 we plot T c as a function 
of t z for the band filling of n = 1.08. As can be seen, the 
three dimensionality suppresses T c , which is in agreement 
at least qualitatively with the experimental finding that 
the intercalation of THF molecules between the layers 
enhances T C J£ 

In fact, the suppression of T c with increasing three di- 
mensionality is a general trend for spin fluctuation me- 
diated pairing as has been studied in ref . 30 ' 31 . Namely, 
since the pairing interaction is large around a certain 
wave vector in spin fluctuation mediated pairing, the 
fraction of the volume where the pairing interaction is 
large in the Brillouin zone becomes smaller as the three 
dimensionality increases. 
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FIG. 7: (color online) Upper two panels : the contour plots 
of the d + id' gap for n = 1.06 (top) and n = 1.16 (middle). 
The model for /3-MNC1 is used. The yellow (light) solid line 
(closed around the Brillouin zone edge) represents the Fermi 
surface. Lower panel : the normalized d + id' gap along the 
Fermi surface. 8 is defined in the top panel. 



V. CONCLUSION AND FUTURE PROBLEMS 

In the present study, we have studied two kinds of 
model for the the doped /3-MNC1. We have found that 
the single band model does not suffice as a minimal ef- 
fective model of the material. On the other hand, from 
the study on the two band model, we have proposed a 
possibility of spin fluctuation mediated d + id'-wave su- 
perconductivity in the doped /3-MNC1. We have shown 
that the relatively high T c is obtained as a characteristic 
feature of the spin fluctuation mediated superconductiv- 
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FIG. 8: T c plotted as a function of t z for the model shown 
in the inset, t = 1.2eV, t' /t = 0.35, A/t = 2.7, U/t = 6, 
n = 1.08. 



ity in a doped band insulator on the HC lattice. We have 
also found that the gap anisotropy on the Fermi surface 
strongly increases upon increasing the doping concentra- 
tion, and the inter layer hopping suppresses superconduc- 
tivity. These results are in qualitative agreement with the 
experimental findings. 

One remaining issue is the curious doping dependence 
of the superconductivity. For Li^ZrNCl, T c shows an 
increase upon lowering the carrier concentration until a 



sudden disappearance of the T c and a superconductor- 
insulator transition is observed^ 5 Recently, this increase 
of T c has been shown to be correlated with the increase of 
the uniform spin susceptibility, which can be considered 
as a support for the present spin fluctuation scenario^ 
In fact, although the spin fluctuation is nearly tempera- 
ture independent in our model (so that the development 
of the spin fluctuation upon lowering the temperature 
need not be observed experimentally) , it is strongly dop- 
ing dependent and is enhanced upon lowering the doping 
rate.— Nonetheless, the origin of the sudden drop of T c 
at a certain doping rate can not be explained within the 
present two band model (see Figj4j) and remains as an 
open problem. Also, the almost constant behavior of T c 
in the heavily doped regime^ is puzzling since for such a 
large doping, the Fermi level hits the second band and a 
large difference in the DOS should be detected (which is 
not the case experimentally). A rigid band picture may 
not be valid in the heavily doped regime, as suggested in 
some experiments^ and band calculations^ 
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